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Abstract 

Gravitational collapse of a spherically symmetric cloud has been extensively studied to inves¬ 
tigate the nature of resulting singularity. However, there has been considerable debate about the 
admissibility of certain initial density distributions. Using the Newtonian limit of the equations 
governing collapse of a fluid with an equation of state p = p{p) it is shown that the density dis¬ 
tribution has to be even function of r in a spherically symmetric situation provided dp/dp ^ 0. 
Implications of this result on formation of strong naked singularities are examined. 


Gravitational collapse of spherically symmetric dust cloud has been extensively studied [1-8] to under¬ 
stand the nature of resulting singularity. Using an initial density distribution which is an even function of 
radial distance, r Newman [3] has shown that the singularity is naked, but weak, when it just forms. By 
dropping the assumption about evenness of density it turns out that with suitable choice of initial density 
and velocity the naked singularity could be strong [4,5]. However, Unnikrishnan [6,7] has argued that such 
a density distribution is not physical when a realistic equation of state with pressure is used. He claims that 
density must be an even function of r. On the other hand, using the special case of polytropic equation of 
state in the Newtonian limit, Jhingan et al. [8] have attempted to show that density is not restricted to be 
an even function of r. Unfortunately, they have not included the continuity equation in their analysis and 
hence their conclusions are not correct. In this work, we show that using equations of fluid dynamics in the 
Newtonian limit with an equation of state of the form p = p{p), (where p is the pressure and p the density) 
density has to be an even function of r in a spherically symmetric situation provided dp/dp ^ 0. 

We consider spherically symmetric collapse of fluid with equation of state p = p{p), which is assumed 
to be continuous and differentiable to any required order. The usual equations of fluid dynamics in the 
Newtonian limit, i.e., the continuity equation and the momentum equation in the Eulerian form are: 
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where v is the radial component of velocity, G is the gravitational constant and 


m= Airr^p dr (3) 

Jo 

is the total mass contained within a sphere of radius r. It may be noted that in [8] only the second equation 
was considered in the Lagrangian form which is equivalent to Eq.2 of [8], while the continuity equation, which 
is the Eq.3 in [8] was neglected. Glearly, the continuity equation is automatically satisfied when considering 
equilibrium situation (v = 0), but while generalizing the Lane-Emden equation to calculate collapse solutions, 
the continuity equation will need to be added to determine the additional unknown, namely, velocity. 

Gonsidering a situation before the singularity is formed, we assume that the solution is regular near the 
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center and expand p, v, p and m in power series about r = 0: 
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Here, m and p are not independent functions but are completely determined by p. Thus we have nii = 
47 rpi /(3 + i) and pi = {dp/dp)pi, etc. There is no constant term in expansion for velocity as v must vanish 
at the center in a spherically symmetric system. 

Substituting these expansions in Eqs. (1,2) and collecting all terms with same powers of r, we can get 
the required equations to determine the coefficients of expansion. Thus considering the terms of in Eq. (2) 
we get 

Pi = ^Pi = 0 (5) 

dp 

Hence, if the equation of state is such that dp/dp ^ 0, then pi = 0 and mi = 0. This condition in physical 
terms simply states that the pressure gradient which is the force, must vanish at the center in a spherically 
symmetric situation. Now, considering terms of in Eq. (1) and using the fact that p\ = 0, we get v\ = 0. 
Similarly, considering terms of in Eq. (2) and using pi = 0, mi = 0 and vi = 0, we get pz = {dp/dp)p 3 = 0, 
and hence p 3 = 0 and m 3 = 0. Again with all these conditions the terms of in Eq. (1) will give V 3 = 0. 
Thus the cubic terms in all expansions also vanish. Continuing this process it is clear that all the odd terms 
in the expansions (4) will vanish. 

Thus it is clear that if dp/dp yf 0, the density is an even function of r, while the radial component of 
velocity is an odd function of r. Therefore, Unnikrishnan’s [6,7] claim that solutions are analytic appears to 
be correct, at least in the Newtonian limit. These results are likely to be valid even for relativistic equations. 
However, it is clear that these restrictions do not apply for the special case of dust {p = 0) considered in 
most calculations of gravitational collapse. Nevertheless, if we think of dust as a limiting case of fluid where 
the pressure gradient is small as compared to the gravitational force, then the density should be considered 
as an even function [7]. 

If we restrict to admissible solutions with p and u/r as even functions then it turns out that the resulting 
solution cannot yield a strong naked singularity [5]. Strong naked singularity can form only when p or v/r 
have some odd terms and as such it appears that these do not arise out of physically admissible initial 
conditions. Even if we restrict to pure dust with no pressure, it turns out that the solutions leading to strong 
naked singularity are unstable (or non-generic) [8,9] and thus physically inadmissible. Although Antia [9] 
considered only the special case of marginally bound system to demonstrate the instability of these solutions, 
the proof is equally valid for a general situation [8]. In the general case, for stability it is essential that the 
solution is stable with respect to small perturbations in both density and velocity. Considering the general 
case of gravitational collapse of spherically symmetric dust cloud, Singh and Joshi [5] have shown that strong 
naked singularity can form when certain quantities Qi and Q 2 defined by them vanish. These restriction 
clearly makes the solutions unstable (or non-generic) since a general perturbation to initial conditions will 
not satisfy these restrictions. Thus even though for any given initial density distribution it may be possible to 
choose initial velocity such that collapse leads to a strong curvature naked singularity [ 10 ], all such solutions 
are non-generic and are not likely to arise in a physically realistic situation. 

In fact the instability of solutions with strong naked singularities is clear from results of Singh and 
Joshi [5], as on p 567 they write “Also, the assumption /2 7 ^ 0 is a natural one, because if the cloud is not 
marginally bound, it will require a fine tuning of the velocity profile with the density profile to make /2 
zero”. In fact, this is precisely the concept of stability or genericity. However, on the next page they consider 
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cases where /a = 0, /4 = 0 etc. If /2 cannot be zero then there is no reason to assume that /a or fi can be 
zero, because that too will require fine tuning. Hence, all these cases, which includes all solutions leading 
to strong naked singularities should not be considered at all. This would confirm the claim of Unnikrishnan 
[6,7] that the naked singularities in physical circumstances can only be weak. 

It may be noted that in all these calculations the singularity is naked only when it just forms. Irrespective 
of initial conditions, gravitational collapse ultimately results in a formation of event horizon and hence black 
holes appear to be the only stable end product of such a collapse. 

We would also like to point out that such non-generic solutions leading to strong naked singularity are 
not likely to be realized in numerical solution [e.g., 11] of equations of gravitational collapse. Since numerical 
solutions are inherently approximate, even if we start with initial values which satisfy the required conditions 
for formation of strong naked singularity, because of errors introduced in numerical solutions these conditions 
will be violated at subsequent time and the resulting solution will ultimately lead to weak singularity which 
results from a generic initial conditions. Similarly, the example given in [8] to illustrate that absence of the 
apparent horizon until singularity formation does not imply that the singularity is naked, is not particularly 
relevant for numerical computations. In generic solutions which are the only ones that would be expected 
in numerical work such behavior is not seen. 

Finally, we would like to point out that even for the case of pure dust (i.e., p = 0), if the density is 
assumed to be an even function of r, then from the continuity equation it is clear that the velocity will be 
an odd function of r (i.e., V 2 n+i = 0, for n = 0,1, 2,...) and once again it is not possible to find any solution 
with strong naked singularity. Thus the solutions with strong naked singularity mentioned in [11] must have 
odd terms in density at a general time t even though at t = 0 the density may be an even function of r. This 
can be seen by writing the coefficient of r in Eq.(2): 

^ -f 4(po'Ci + Pi^^o) = 0. (6) 

Thus if at t = 0, Pi = 0 but vi ^ 0, then dpi/dt ^ 0 and hence the density will pick up odd terms. Such 
solutions have no significance, as the density is an even function of r at t = 0 but does not remain so at any 
other time. Thus the choice of time t = 0 at which the collapse is supposed to start has to be very special 
to yield such solutions. In general, it should not matter which epoch is chosen as t = 0 and hence, if there 
is any reason for the density to be an even function of r it should be so at all times t, in which case velocity 
must be an odd function of r and it is not possible to find any solution with strong naked singularity. 
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